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Abstract 

A model for the adsorption of a binary mixture on a one-dimensional infinite lattice with nearest 
neighbour cooperative effects is considered. The particles of the two species are both monomers but 
differ in the repulsive interaction experienced by them when trying to adsorb. An exact expression 
for the coverage of the lattice is derived. In the jamming limit, it is a monotonic function of the 
ratio between the attempt frequencies of the two species, varying between the values corresponding 
to each of the two single species. This is in contrast with the results obtained in other models for 
the adsorption of particles of different sizes. The structure of the jamming state is also investigated. 
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I. INTRODUCTION 



The adsorption of particles on a solid substrate is a process occurring in many natural 
phenomena and that, consequently, has been extensively studied both experimentally and 
theoretically The simplest kind of models for adsorption are those defined in terms 
of random sequential adsorption (RSA) processes, in which particles of a given kind are 
adsorbed at random positions on the substrate. When the adsorption rates depend on the 
state of the neighborhood of the position at which they take place, the process is known 
as cooperative sequential adsorption (CSA). A lattice version of the CSA process is the 
monomer filling with nearest-neighbor (NN) cooperative effects ^,|2|- These lattice models 
have been extensively used to study the monolayer growth of monodisperse particles. In 
particular, exact analytical solutions have been derived for the one-dimensional case 

The deposition of a mixture of different kind of particles has also received some attention 
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, but the progress made in its understanding is small as compared 
with the one-component problem. Here, a one-dimensional lattice model for a binary mixture 
of monomers with NN cooperative effects is considered and its analytical solution is found. 
Besides their theoretical interest, one-dimensional models are also useful for the description 
of some polymer reactions Q]. A controversial point in the adsorption of a mixture is whether 
it covers the substrate more or less efficiently than either of the species separately. Some 
results seemed to indicated than the former was true for lattice models 0,0], while the latter 
applied for continuum models Nevertheless, a (continuum) random car parking model 



has been found to be consistent with the lattice models in this respect jlO, For the 
lattice model considered here, it will be shown than the jamming coverage of the mixture is 
always between the coverage associated with each of the two components. This appears to 
be a direct consequence of the NN cooperative effects taken into account in the model. 

It is interesting to note that lattice models with NN cooperative effects have also been 
considered in the context of Ising models. They were introduced by Fredrickson and Ander- 
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14j to study structural relaxation in glassy systems. In these models, a given 



spin can flip only if its nearest neighbors are in a certain subset of all their possible config- 
urations and, because of this, they were termed "facilitated" Ising models. In particular, in 
the so-called nSFM, one spin can flip only if at least n of its NN spins are orientated against 
the external field. It is easily verified that the zero temperature limit of these models is the 
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spin representation of the hole-particle (monomer) models mentioned above, i.e. both are 
related by a single change of variables. Interestingly, the one-dimensional 1SFM has also 
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been used as a model to describe compaction proceses in vibrated granulates 

In the model considered in this paper, the two kind of particles composing the mixture 
differ in their dynamics. Particles of one of the species obey RSA with NN exclusion or 
blocking i.e. they can occupy an empty site with the constraint that both of its NN 
must be empty. On the other hand, adsorption of one particle of the other species on an 
empty site only requires that one of the NN be empty. Then, the dynamics of both kind of 
particles is cooperative, but the "facilitation" rules are different. A possible interpretation 
of this form of competitive adsorption is that particles are of the same size, but the effective 
repulsive interactions coming from the adsorbed particles are much stronger for one species 
than for the other. Then, the repulsive effect depends on the kind of particle trying to 
adsorb and on whether the NN sites are occupied or empty, but not on the type of particle 
actually occupying them. 

In the next Section, the model will be formulated in detail and the master equation 
describing its dynamics will be written down. The dynamics of the adsorption processes 
is studied in Sec. IHIl where explicit analytical expressions for the time evolution of the 
densities of each component are derived. For very large times, the system gets stuck in 
a jammed configuration, where no more adsorption events are possible. The asymptotic 
values of the densities of both species, as well as the total coverage, are also calculated as a 
function of the adsorption rate constants of the two species. It is found that the maximum 
asymptotic total coverage is always smaller than the one corresponding to the least repulsive 
particles. In order to get a deeper understanding of the stationary, jammed, state reached 
by the lattice in the long time limit, some correlation functions are analysed in Sec. II VI The 
jammed states are characterised by having all the holes isolated, i.e., surrounded by two 
particles. Therefore, we will study the structure of the lattice around these isolated holes, 
analysing the probability of finding them between all the possible configurations of their NN 
sites. Section contains a short summary of the main results and conclusions of the paper. 
Finally, some calculations are presented in an Appendix. 
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II. THE MODEL 



We consider two kinds of particles, A and B, that can be randomly adsorbed on a one- 
dimensional lattice. The adsorption processes on an empty site are restricted by the config- 
uration of its NN, as follows: 

(i) Particles A can be adsorbed on an empty site only if at least one of its NN is also 
empty. 

(i) Adsorption of a particle B on a site requires that both of the NN are empty. 

Let us introduce two occupation numbers rii and for each site i of the lattice. If there 
is a particle A at site i, it is rii = 0, otherwise rii — 1. Similarly, mj = 1 if there is a particle 
B at site i, and rrii = if there is not. Of course, since there can be at the most one particle 
at each site, it is (1 — 7ij)(l — = for all i. A lattice configuration is specified by the 
sets of occupation numbers of the two species {n, m}, where n = {rii} and m = {rrii}. 

The dynamics of the system is defined as a Markov process, so it is enough to specify 
the master equation for the probability distribution of the configurations, p(n, m, t). As 
said above, the elementary processes are the adsorption of particles A and B on an empty 
site i. In these events, the configuration of the system changes from (n, m) to (RiU, m) or 
to (n, Rim), depending on whether the adsorbed particle corresponds to species A or B. 
Here Ri is the operator changing the occupation number rii or m i into 1 — rii or 1 — n^i, 
respectively. Thus, the master equation reads 

^-p(n,m,i) = [Wj(n, m\Rjn, m)p(Rjn, rn. t) — Wj(Rjn, m\n, m)p(n, m. t)\ 

i 

+ ^ [Wj(n, m\n, Rim)p(n, Rim, t) — Wi(n, Rim\n, m)p(n, m, t)] . (1) 

i 

where 

OL 

Wi{Rin, m\n, m) = -riivrii (nj_imj_i + n i+1 m i+1 ) , (2) 

Wi(n, Rim\n, m) = Prii-imi^riimirii+imi+x . (3) 

The above transitions rates describe the adsorption of a particle A and B on site i, respec- 
tively, consistently with the rules (i) and (ii). The constant parameters a and (5 characterise 
the attempt frequency for each kind of adsorption process. 

If a = or (3 = 0, the lattice will be occupied by particles of only one kind, and the 
dynamics of the model reduces to the one-component monomer filling with NN cooperative 
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effects, whose solution is known and discussed in detail by Evans in Ref. When both 
rates are nonzero, we have the competitive adsorption of the two species: particles B need 
both of the NN sites of a given hole to be also empty in order to adsorb, while one e mpty 
NN site is enough for a particle A. In the corresponding facilitated Ising picture 12l ll3L ll4|. 



if (3 = we have the one- dimensional 1SFM, while for a = the one-dimensional 2SFM is 
recovered. 

We will restrict ourselves to translationally invariant states. This requires the consid- 
eration of consistent boundary conditions, e.g. periodic ones. In this case, the density of 
particles A at a given time t is given by 

p A (t) = i-(m) t , (4) 

where the angular brackets denote average with the probability distribution p(n, m, t) and 
the right hand site does not the depend on the site i chosen. Similarly, the density of B 
particles is 

p B (t) = 1 - (mi) t . (5) 
The total coverage of the line 9(t) is then 

9(t) = p A (t) + p B (t) . (6) 

It verifies < 9{t) < 1, since it is impossible to fully fill the line, as it will be discussed 
in detail latter on. A particularly interesting property is the asymptotic coverage 9j at 
jamming, i.e., once the system gets stuck in a configuration in which no more adsorption 
events are possible. In the model considered here, this happens when all the holes are 
isolated. Formally, we can write 

9 j = lim 9{t) . (7) 

t— >oo 

It is useful to consider the set of moments 

F r (t) = (nimin i+1 m i+1 . . . n i+r m i+r ) t , (8) 

for all r > 0. They give the probability of finding r + 1 consecutive holes or, equivalently, 
the density of clusters of at least r + 1 empty sites. In particular, F (t) gives the density of 
holes and, therefore, it must be 

F (t) + p A (t) + p B (t) = 1. (9) 



The above equation can also be written as 

((1-70(1-^ = 0, (10) 

expressing that a site can not be simultaneously occupied by a particle A and a particle B. 
Besides, between each pair of particles B there must be at least one empty site, due to the 
adsorption rule for particles A. Therefore, 

F Q (t) = 1 - p A (t) - p B (t) >p B (t). (11) 



III. ANALYTICAL SOLUTION OF THE DYNAMICS 

Here, the consequences of the dynamics defined in the previous section will be analysed. 
First, the time evolution of the densities of both species will be investigated. Equations for 
them are readily obtained from the master equation (JTJ), 

-p A (t) = aF 1 {t) , (12a) 

j t p B {t) = PF 2 (t) , (12b) 

where Fi(t) and F 2 (t) are the densities of clusters with at least two and three consecutive 
holes, respectively, defined in Eq. (jHJ). Also, a hierarchy of equations for all the marginal 
probabilities F r (t) is derived from the master equation, 

j t F (t) = -aF x {t) - /3F 2 (t) , (13a) 

jF r {t) = -(a + 2(3)F r+1 (t) - [ar + (3{r - 1)] F r (t) , r > 1 . (13b) 

Combination of equations (|T2"|) and ()13aj) implies that F (t) + pAif) + PBif) is in fact an 
integral of motion, as required by Eq. ©. In order to solve the hierarchy of equations 
()13b|) . it is convenient to introduce the generating function 

oo r 
X 



G(x,t) = J2^F r+1 (t) J (14) 

such that 



T 
r=0 



The generating function satisfies a linear first-order partial differential equation, namely 

d t G(x, t) + [a + 2(3 + (a + f3)x] d x G(x, t) = -aG(x, t) , (16) 



which has to be solved with the initial condition 

oo r 

G(x,0) = G (x) = £^F r+1 (0) . 



r=0 



By using standard techniques, it is obtained that 

a + 2(3 



G(x,t) = G c 



+ 



-at 



a + (3 \ a + (3 

Then, from Eq. (JT5j) we derive the whole set of moments F r (t) for r > 1, 



F r {t) = G, 



(r-l) 



a + (3 

where we have introduced the notation 



a + 2(3 , (a + 2(3 
a + (3 



+ 



-{a+0)t 



e -art e -f3(r-l)t 



G 



x 



d r G {x) 
dx r 



(17) 



(18) 



(19) 



(20) 



In the long time limit t — > oo, all the moments -F r (t) such that r > 1 vanish, provided that 
a^O. This is easily understood: if adsorption of particles A is possible, the system evolves 
until it gets stuck in a metastable configuration such that all the holes are isolated, i.e., 
there are no pairs of consecutive holes. On the other hand, if a = 0, so that only adsorption 
of particles B is allowed, i^(oo) ^ while F r (oo) = for r > 2. In the long time limit 
those configurations with at most two consecutive empty sites are jammed since, in order 
to adsorb on a given site, particles B need both of the nearest neighbours being empty. 
Therefore, the limit a — > must be handled with care in the present model, as the jammed 
configurations are different for a = and for a — > but not identically null. 

The evolution equations for the densities of particles pa($) an d Pb(£) are obtained by 
substituting Eq. (JTHJ) into Eqs. ([T2|). 



d , . 



d 
It 



PB(t) 



a Gq 



a + 2(3 a + 2(3 



+ 



a + (3 a + (3 
a + 2(3 a + 2(3 



+ 



-(a+/3)t 



-<\t 



-{2a+f3)t 



(21a) 
(21b) 



a + (3 a + (3 

These equations must be integrated to get explicit expressions for the densities pA,(t) and 
Pb{£)- In order to do this, we have to specify the initial condition for the generating function 
Go{x) or, equivalently, the complete set of marginal probabilities F r (0). 

For the sake of concreteness, we will consider that the lattice is empty at t — 0. This is 
the usual initial state in adsorption studies. Therefore, 



^(0) 



Vr > 0, 



(22) 



and, as a consequence, 

From Eq. (JUJ) we get 

F r (t) = exp 



r=0 



a + (3 a + (3 
for r > 1. In the long time limit t — > oo, 



a + 2/5 a + 2/5 r„,flw 

+ _£ e -(«+p)* _ art _ _ 1)£ 



F r (oo) = lim F r (t) = , Vr > 1 



(23) 



(24) 



(25) 



provided that a^O, consistently with the above general discussion. As already indicated, 
for a = the adsorption of particles A is impossible, and, in general, -Fi(oo) ^ 0. Use of 
Eq. flU into Eq. fl2U) yields 

a + 2(3 , a + 2/5 _ (a+/3)t 



eft 



Pa(^) = a exp 
p s (i) = /5exp 



, H 7T e 

a + (3 a + (3 

a + 2(3 a + 2(3 _ 



at 



+ 



-(a+/3)t _ 



a + (3 a + (3 
These equations can be integrated, with the result 

1 l + 2r 

-e 1 + r 



a + 2(3)t 



PA(t) 

PB(t) 



1 + r 



r l + 2r 



-(<*+|3)t 



f l + 2r 

e !+ r 

1 + r j ( a+ p) t 



_J_ l + 2r 

cm it !+»■ e 1 + r 



where we have introduced the ratio of the attempt frequencies of adsorption 

a 



(26a) 
(26b) 

(27a) 
(27b) 

(28) 



Integration by parts of Eq. (|27b|) gives an explicit relationship between both densities pAif) 
and p B (t), 

P^-TT^-i 1 '^' 6 ^ 1 ''''"''^-^}- (29) 

Let us analyse the long time limit of the densities, when the system gets jammed in 
a metastable configuration and no more adsorption events are possible. In this limit, the 
densities tend to the values 



p A {oo) 
Pb(oo) 



1 l + 2r / r_ l + 2r 

-e 1 + r / duu !+ r e 1+r 



1 + r 



1 + r 



r l + 2r / _J_ l + 2r 

-e *+r / auu 1 + r e 1 + r 



(30a) 
(30b) 
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which are related by 



T 



as long as a 7^ 0. In the long time limit, the tendency to these asymptotic limits are 

Pa{oo) - PA {t) ~ e-^e- at , (32a) 
Pb(oo) - PB (t) ~ -^e-^e-^ 4 . (32b) 

Both behaviours are exponential, but the characteristic time (2a + /3) _1 for particles B is 
smaller than the characteristic time a -1 for particles A. Thus, the total coverage of the line 
0{t), defined by Eq. (0), verifies 

9j - 9(t) = PA (oc) - PA (t) + p B (oo) - p B (t) = e-^e- at + O [ e ~^ + ^] , (33) 

where 9j is the total coverage at jamming, given by Eq. (JJJ). The asymptotics of the 
coverage is dominated by the contribution of particles A. This is easily understood, since as 
time increases a configuration will be reached where no more adsorption of particles B are 
possible. On the other hand, particles A can still be adsorbed on the lattice. Our results 
show that this situation will be found for times or the order of (2a+f3)~ 1 , for which ps{t) has 
already decayed to ps(oo) while p A {t) is still evolving. This exponential approach near the 
jamming limit is characteristic of any adsorption process in a one- dimensional lattice. Power 
aw convergence, like Feder's t" 1 law arises when continuous deposition is considered 

mm 

In the limit r — > 00 (a <C /?), Eq. (j30b|) leads to 

limp B (oo) = i(l-e- 2 ) , (34) 

and, therefore, 

lim pa(oo) = 1 — 2 lim p B (oo) = e~ 2 , (35) 

1 — »oo r— >oo 

which is nonzero. On the other hand, for a = we know that p A (t) = for all t. This shows 
again the singularity of the case a = 0. For any a ^ 0, the system gets stuck in configu- 
rations characterised by having all the holes isolated, while for a = those configurations 
having, at most, two consecutive holes are also metastable. The adsorption of particles A 
on a previously jammed configuration of particles B is analysed in the Appendix. In that 
situation, the asymptotic density of particles A, Eq. (|Allj) . equals Eq. (J33J). For a — > but 
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FIG. 1: Plot of the asymptotic densities pa(oo) (crosses) and pb{oo) (triangles) as a function of 
the adsorption rates ratio r. Note that lim r -_ >0 c Pa(oo) / 0, as expressed by Eq. (|H5jl. 

nonzero, the system first reaches, over a time scale of the order of f3~ , a jammed configura- 
tion with only particles B. Afterwards, for times of the order of a -1 ^> particles A are 
adsorbed on this state, leading to an asymptotic density of particles A given by Eq. (J35J). 
The asymptotic total coverage, defined in Eq. (JJ|), in the limit case we are considering is 

lim 9j = lim \p A {oo) + pb(oo)] = - (l + e" 2 ) , (36) 

r— >oo r— >oo ^ 

which also equals Eq. (|A12J) . as expected on the basis of the discussion above. 

The limit r — > 0, i.e., /3 — > does not present any kind of singularity. For r = is is 
readily obtained that 

limp^oo) = 1 — e _1 , limps(oo)=0. (37) 

The system only contains particles A, and the result agrees with the one obtained for the 
one-component system 

As it would be expected on physical grounds, the asymptotic density of particles ps(oo) 
is a monotonic increasing function of r, while p^(oo) decreases with r. A plot of both 
asymptotic densities as a function of the adsorption rates ratio r is shown in Fig. ^ The 
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jamming coverage can be written in terms of p^oo) by using Eq. (|31jh 

6j = PA {oo) + PB {oo) = y~^ + TT^ pA (°°^ ( 38 ) 

which is a decreasing function of r, varying from 

limflj = 1 -e- 1 ~ 0.63 (39) 

to 

lim flj = - (1 + e~ 2 ) ~ 0.57. (40) 

Figure shows the total coverage 8j as a function of r. The maximum value of the coverage 
is obtained for r — > (/3 — > 0), i.e., when the adsorption of particles B is forbidden, and 
there are only particles A in the system. This is reasonable on physical grounds: due to the 
facilitation rule, the domain between two particles B cannot be completely full of particles A. 
Therefore, the density of particles B is a lower bound for the density of holes. As the density 
of particles B increases with r = f3/a, the asymptotic density of holes F (oc) also increases 
with r and the jamming coverage 6j decreases, since 6j = 1 — -Fo(oo). This behaviour is in 
contrast to all the previous results for RSA in continuum parking models and in 

lattice models In these works, it is reported that the substrate is covered either less 

or more efficiently by a binary mixture than by a single species. The model presented here, 
however, leads to binary coverage lying between that of the two monodisperse cases. This is 
a consequence of the cooperativity of the adsorption processes taking place in our system. 



IV. CORRELATIONS 



In order to get a more complete description of the steady state, we have also investigated 
some correlation functions. In the steady state, holes are isolated, i.e., surrounded by two 
particles. In particular, we will be interested in the structure of the lattice surrounding these 
isolated holes. In the previous Section, the time evolution of the total hole density -Fo(^) 
has been found, and here the density of holes for a given configuration of their two nearest 
neighbours will be studied. In the asymptotic stationary state, they can be both particles 
A (AA), one particle A and one particle B (AB or BA), and two particles B (BB). Let us 
define the density 

®o A {t) = ((! - ni_i)Trii_ini77ii(l - n i+1 )m i+l ) t , (41) 
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FIG. 2: Plot of the asymptotic total coverage 6j as a function of r. 

that corresponds to the density of empty sites surrounded by two particles A. In a similar 
way, we introduce 



$o S (0 = ((! - ni-^mi-inimirii+iil - m i+1 )) t , 



(42) 



$o = - m i _i)7i i 7n i (l - n m )m m ) t , 



(43) 



$0^(0 = - mi_i)njmin i+ i(l - m m )) t , (44) 

corresponding to the density of holes between a particle A on the left and a particle B on 
the right (§q B ), etc. In the jammed state, 

F (oo) = $r (oo) + $^(oo) + < A (oo) + < B (oo) . (45) 

Nevertheless, an analogous expression does not hold for arbitrary time t, since then it is 
possible to find configurations with two adjacent holes. This is clear from the fact that the 
moments F r (t), with r > 1, do not vanish except in the limit t — > oo. Then, the asymptotic 
fraction of holes between two particles A at jamming is 



x 



AA 



$0^(00) 

F (oo) 



(46) 



12 



Similarly, we define 



= (47) 



as the relative density of holes between a particle A to the left and a particle B to its right, 

^ = T7=P (48) 

F (oo) 

which gives the relative density of holes between a particle B to the left and a particle A to 
its right, and 

*?" = T7=r- (49) 

F (oo) 

for the fraction of holes between two particles B. 

The evolution equations for these moments are again readily obtained from the master 
equation. For the density of holes between two particles A, one gets 

j t ^\t) = a [*}(t) + *i{t)] , (50) 

where 

§f(t) = (nimini +1 m i+1 . . . % r _im i+r _i(l - n i+r )m i+r ) t , (51) 

i.e., it is the probability of finding r consecutive holes from site % onwards, and a particle A 
on site i + r. These probabilities obey the hierarchy of equations 



V = - (f + P) <i - Hr - 1) + P(r - 2)] + | (F r + F r+1 ) , (52) 



dt 



for r > 2. For an arbitrary initial condition, this hierarchy of equations can be solved by a 
generating function method analogous to the one used to integrate the evolution equations 
for F r (t), in spite of the inhomogeneous term proportional to F r + F r+ x. Nonetheless, for 
the initially empty lattice configuration, it is simpler to assume that 

^(t)=U(t)e~ [a{r ~ 1)+P{r - 2)]t ■ (53) 

This time dependence is suggested by the own hierarchy of equations (|52|) and, also, by the 
time dependence of the inhomogeneous term, which is precisely of this form, as follows from 
Eq. ((2H) . By substituting Eq. (jKHjl into Eq. (JE2J), a differential equation for £a(£) follows, 



dt 



(| + (3) e-^U + | exp + ^§ e ~ (a+P)t e ~ {a+P)t I 1 + e ~ (a+/3) *] 



a + (3 a + P 

(54) 
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This is easily integrated, with the result 



2a(5 



[a + 2(5)- 



a 



a + 2(3 



exp 



cxp 



a + 2(5 a + 2(5 



+ 



2(a + (5) 2(a + (5) 
a + 2p + a + 2(5 c _ {a+0)t 



a 



a + 2(5 



_ e -(a+p)t 



(55) 



a + (5 a + (5 

Thus, we have identified all the moments 3>Jr (0 and, by using Eq. 1)50)1 . obtain the following 
closed evolution equation for <3>q 



d 



^\ t )=aUit)e- at [l + e-^] . 
Since the lattice was initially empty, § A A {t = 0) = 0, and 



(56) 



2r 



l + 2r 



;i +r)(l + 2r) 2 
1 



l + 2r 



r 1 ■ : / duw i+-(l + w)e 2 ( 1 +''>" 

(t»+|8)« 

1 



+ 



u eW u , (57) 



l + 2r 

l + r)(l + 2r) y e -(a+«* V ; Vl + 2r 

where r was defined in Eq. (128)1 . Therefore, the asymptotic density of holes surrounded by 



two particles A is 



2r 



(1 + r)(l + 2r) 2 
1 



_ l + 2r / l+2r 

e 2(T+7T / — (l + u ) e 5TL+?7 u 



, e !+'- / duw !+'• (1 + u, , 
l + r)(l + 2r) J V '\l + 2r 



\ \-+2r 

u )e^ u . (58) 



In the limit r — > 0, this equation reduces to 



lim$^ A (oo) = e- 1 [ du (1 + u)(l - w)e" = e" 1 = lim F ( 



OO . 



(59) 



This is readily understood, as the limit r —>■ corresponds to (3 — 0, in which no adsorption 
events for particles B are possible. On the other hand, in the limit r — > oo, 



lim $o A (°°) = °- 



(60) 



This result may appear as nontrivial, since it must be stressed that when r — > oo there 
is a finite, nonvanishing, density of particles A in the steady state, given by Eq. ([35)1. 
Nevertheless, as discussed below that equation, the adsorption of particles A takes place for 
very long times t = 0(a^ 1 ) ^> on sites with at most one empty nearest neighbour. 
As the latter must be next to a particle B, it follows that <& AA vanishes for r — > oo. In 
Fig. El the fraction of holes between two particles A x AA , defined in Eq. (jjBj) , is plotted 
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FIG. 3: Plot of the relative density Xj A (oo) (squares) as a function of the adsorption rates ratio 
r (in a logarithmic scale). 

as a function of r. It decreases monotonically from unity for r = to zero for r — > oo, as 
expected. 

Let us now analyse the density of holes between two particles B, $^ B (t). From the master 
equation, it is obtained 



d 
dt 



BB 



2/3$ 



B 
3 ) 



(61) 



where now 

$f (t) = (nimin i+1 m i+1 . . . n i+r _im i+r _in i+r .(l - m i+r )) t , (62) 

i.e., it is the probability of finding r consecutive holes from site i onwards, and a particle B 
on site % + r. Note the analogy of the notation used here with the one introduced for the 
moments <3>^ defined in Eq. (151 J) . These moments obey the hierarchy 

j$r = - (f + P) *£u - [a(r - 1) + 0(r - 2)] $f + /3F r+1 , (63) 
for r > 2. Again, for the initially empty lattice we are considering, the family of solutions 



-[a(r-l)+P(r-2)]t 



(64) 
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can be considered. Substitution of the above expression into Eq. yields a closed differ- 
ential equation for namely 



Jt iB 



(a+f3)t 



£ B + (3 exp 



a + (3 a + (3 



-2(a+P)t 



(65) 



which can be integrated by standard techniques, obtaining 



4/?(q + g) 



(a + 2/?) 2 
2a/? 
"(a + 2/3) : 



a + 2/3 | a + 2^ c _ (a+/3)t 



exp 



a + (3 a + (3 
® + W + a + 2(3 c _ (a+0)t 



a + 2(3 
2(a + (3)' 



-(a+/3)t 



(66) 



2(a + (3) 2(a + (3) 

This provides the expressions for all the densities $^ (r > 2), from Eq. (J64j) . Taking into 
account Eq. (JBT)) and the initial condition, it is found 



8r 2 



l + 2r) : 



l + 2r 

-e 1 + r 



4r z 



-(a + /3)t 



l + 2r 



P 2(l + r) 

;i + r)(l + 2r) 2 ./e-(»+«« 



1 + 2r l+2r . 

1 — — -w e !+ r 

2(1 + r) 

1 l + 2r 



(67) 



In the long time limit, 



2 /■! 

r _i+2r 
e 



(l + 2r) 2 



duu 1 ^ 



o 



1 + 2r 

1 ; -U 

2(1 + r) 



l + 2r 

e 



4r 2 



(l + r)(l + 2r) 2 7 
If the adsorption of particles B is forbidden, (3 — ► or r — > 0, the expected result, 

lim< B (oo) = 0, 

i — >o 

is obtained. On the other hand, when no adsorption of particles A, a — > or r — > oo, 



lim < s (oc) 

r— >oo 



1 



2e z / du(l - u)e"" - - -< 



2u — _ _ _ f - 2 
" 2 2 



(68) 



(69) 



(70) 



Then, $^ B (oo) does not equal the total density of holes Fo(oo) in this limit. This is due to 
the fact that for a^O, although very small, there is a non vanishing density of particles A, 
given by Eq. (J35j) . The consistent result is 



lim [F (oo) - $^ B (oo)] = lim p A (oo) 



(71) 



As we have already discussed above, in the limit r — > oo particles A are adsorbed once the 
system has reached a jammed configuration in which only particles B have been adsorbed 
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FIG. 4: Plot of the concentration function of r (in a logarithmic scale). Note that Xj 

does not tend to unity for r — » oo. 



BB 



(a = 0). These intermediate configurations are characterised by having clusters of holes 
involving at most two consecutive sites. Therefore, particles A will be adsorbed in one of 
the empty sites of those clusters with two adjacent holes, i.e., 



lim [< B (oo) +p A (oo)] = lim F (oo) 



(72) 



because pa(oo) equals the density of holes between a particle B and another hole in the 
intermediate "metastable" state. In Fig. 01 the fraction of holes between two particles 
B, Xj B , is plotted as a function of the adsorption rates ratio r. It is seen that Xj B is a 



monotonically increasing function of r. It must be noted that lim^oo x BB ^ 1, since in that 
limit F BB (oo) does not equal the total density of holes F (oo), as expressed by Eq. (|72*jl . 
More concretely, 

rUL! '^ ] 1 " 3r " 2 (73) 



lim xj B 

r—*oo 



r-»oo Fq{00) 



Taking into account Eq. (|45j). it is possible to obtain also the density of holes between 
one particle A and one particle B, independently of their relative positions 



(74) 
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FIG. 5: Plot of xj B {oo) as a function of r (in a logarithmic scale). A maximum occurs for r ~ 2. 

The fraction of holes between two distinct particles Xj B = x BA is plotted in Fig. E| It 
has a non-monotonic behaviour, exhibiting a maximum for r ~ 2, i.e., in the region where 
the attempt rates of adsorption a and (3 are of the same order of magnitude. For r — > 0, 
Xj S — > 0, since no particles B are adsorbed on the lattice. On the other hand, in the limit 
r^ooa "minimum" adsorption of particles A is required on those groups of two consecutive 
holes left by the previous adsorption of particles B, as already discussed. Then, 

lim xj B = lim ° / / = lim H *\ \ = — , (75) 

rwoo r^oo Fo(oo) r^oo 2^0(0X3) 1 — 6 2 

since in that limit + ^o A = 2$q B should equal pa, as expressed by Eqs. (J7T|) or (|72*j) . 



V. CONCLUSIONS 



The objective here has been to study a one-dimensional model for the adsorption of a 
binary mixture on a lattice with nearest neighbour cooperative effects. Both species are of 
the same size (monomers), but they differ in the strength of the NN repulsive interactions, 
as measured by the number of NN empty sites required for a particle to adsorb. Emphasis 
has been put on the jammed configuration reached by the system in the long time limit. 
The asymptotic coverage of the lattice by the mixture differs from previous results both 
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'or the random sequential adsorption of mixtures on lattices and continuous substrates 
IflflHEiL The mixture a!wa yS cove, the lattice m o r e efficient,, than the species of 
particles feeling stronger NN repulsion, and less efficiently than the species experiencing a 
weaker one. In fact, the jamming coverage is a monotonic function of the ratio between the 
attempt rates of adsorption for the two species. This is a consequence of the cooperativity 
of the adsorption events considered in our model. 

In the long time limit, the tendency to the asymptotic value of the coverage is dominated 
by the least cooperative species. This is similar to the behaviour observed in models with 
species of different sizes, in which this approach is dominated by the smaller species jj, s[ 
Nevertheless, in our model the approach is exponential and, therefore, Feder's 



8 Jeuitssui < 

[9, 10, [n|. 



t^ 1 power law ^| does not hold. This is a consistent result, since exponential behaviour 
is characteristic of any d-dimensional lattice deposition, evolving to a power law when a 
continuum deposition limit is introduced 

The structure of the jamming configuration has also been investigated by computing some 
spatial correlation functions. Here, the interest has been put on the nature of the holes, i.e., 
the kind of particles surrounding them. The relative densities of holes between two identical 
particles, of either of the two species, show a monotonic behaviour as a function of the 
adsorption rates ratio. On the other hand, the fraction of holes between two particles of 
different species displays a maximum in the region where both attempt rates of adsorption 
are of the same order of magnitude. This may be a relevant information in order to determine 
the properties of the adsorbed layer in a given physical system. 

The model discussed here is closely related to the so-called facilitated Ising models, first 
introduced by Fredrickson and Andersen 

0, 0, Q. The main characteristic of these 

systems is that, although their thermal equilibrium properties are trivial, their dynamics 
is highly cooperative and very slow. A recent review on these kinetically constrained Ising 



models can be found in Ref. [2l| . 
systems as structural glasses l l 



he y ar e appropriate and have been used to model p 



22l I23L |24( or dense granulates 



23,121 



rysical 



111 



which the structural rearrangement of a given region may be slowed down or even blocked 
by the configuration of the surroundings. 



re 



Irreversible cooperative adsorption processes have also been used as models 
axation between two shakes in vibration experiments with granular systems 



or the 
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25 



rcc 



26 



271 . 128L l29f . Let us consider an horizontal section of a real granular binary mixture, near 
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the bottom of its container. During the free evolution of the system, i.e., only under the 
action of gravity, particles can only go down, as long as there is enough empty space in 
their surroundings. The total density of particles in the layer grows until the hard-core 
interaction prevents more movements of particles, and a metastable (mechanically stable) 
configuration is reached. Different kind of grains may need also different amounts of free 
space in their surroundings to be adsorbed on or desorbed from the layer. In our system, 
these kinetic constraints are modeled by the different "facilitation rules" for the adsorption 
of both species. On the other hand, when the system is submitted to vertical vibration, 
particles can go up, decreasing the density in the layer. This vibration dynamics can be 
modeled by allowing desorption events. In this way, a simple model to analyse granular 
segregation phenomena in a binary mixture can be formulated jl7|. 
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APPENDIX A: ADSORPTION OF PARTICLES A ON A JAMMED CONFIGU- 
RATION OF PARTICLES B 

In this Appendix, we will briefly analyse the dynamics of a system in which particles 
A adsorb on a jammed state corresponding to a previous adsorption of particles B. The 
adsorption of B particles is equivalent to the one- dimensional RSA with NN exclusion or 



blocking, whose solution is well-known 
of holes reads 
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. In the jammed state of particles B, the density 



^o(oo) = ^(l + e- 2 ) , (Al) 
and the density of pairs of consecutive holes is 

Fi(oo) = e" 2 . (A2) 

Since the jammed state of particles B is characterised by having at most two consecutive 
holes, 

F r (oo) = for all r > 2. (A3) 
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The asymptotic density of particles B directly follows from the density of holes, 

p B (oo) = 1 - F (oo) = 1 (1 - e- 2 ) . (A4) 

Let us now consider that this jammed state for the system of particles B is the initial 
state for the adsorption of particles A. Therefore, in terms of the moments F r (t) we will 
have 

F r (0) = 0,Vr >2; ^(0)^0, F (0)^0, (A5) 

being Fq(0) and -Fi(O) the initial density of holes and of pairs of consecutive holes, given by 
Eqs. (I A 1(1 and (jA2|) . respectively. Besides, the initial densities of particles are 

Pa(0) = 0, p B (0)^0, (A6) 



where Pb(0) is given by Eq. (|A4J) . 

In order to solve this problem, the formalism developed for the complete system in Sec. 
TTT1 making (3 = 0, can be used. The initial generating function is, 



oo 

X' 



Go(x) = J2^Fr + i(0)=F 1 (0), (A7) 

r=0 

which is independent of x. Then, Eq. (J19)) yields 

F 1 (t) = F 1 (0)e-°*, (A8) 

i.e., the density of pairs of consecutive holes decays exponentially. Of course, F r (t) = for 
all r > 2. The evolution equations for the densities (|2"Tj) are very simple, 

j t p B (t) = , j t p A (t) = aF^)e- at , (A9) 

leading to 

p A (t) = Fi(0) (1 - e- at ) . (A10) 

This result reflects that the initial configuration is a jammed state of the particles B and, 
therefore, the holes are either isolated or in groups of two consecutive holes, whose density is 
F\(0). Particles A can be adsorbed on any of the 2F\{t) sites which are next to a hole, with 
a rate a/2. As all of these processes are independent, the density of hole pairs F±(t) decays 
exponentially as exp(— at), and the density of particles A also increases exponentially until 
it reaches its steady value 

p A (oo) = F 1 (0)=e- 2 , (All) 
21 



with a characteristic time r = a 1 . Note that this is precisely the limit of p^oo) when 
r = (3 /a — > oo, Eq. Besides, the total coverage of the line will be 

6j = p^(oo) + p B (oo) = ~ (1 + e" 2 ) , (A12) 

which agrees with Eq. (j36j) . 
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